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ABSTRACT 

When measuring diameters of partially resolved sources like planetary nebulae, Hii 
regions or galaxies, often a technique called gaussian deconvolution is used. This tech- 
nique yields a gaussian diameter which subsequently has to be multiplied with a con- 
version factor to obtain the true angular diameter of the source. This conversion factor 
is a function of the Full Width at Half Maximum (fwhm) of the beam (in the case of 
radio observations) or the point spread function (in the case of optical observations). 
It also depends on the intrinsic surface brightness distribution of the source. 

In this paper the theory behind this technique will be studied. This study will be 
restricted to circularly symmetric geometries and beams. First an implicit equation will 
be derived, from which the conversion factor for a given surface brightness distribution 
and beam size can be solved. Explicit expressions for the conversion factor will be 
derived from this equation which are valid in cases where the beam size is larger 
than the intrinsic size of the source. A more detailed discussion will be given for 
two simple geometries: a circular constant surface brightness disk and a spherical 
constant emissivity shell with arbitrary inner radius. The theory is subsequently used 
to construct a new technique for determining the fwhm of an arbitrary observed 
surface brightness distribution. 

Usually the fwhm of the source and beam are measured using gaussian fits, but 
second moments can also be used. The alternative use of second moments in this 
context is studied here for the first time and it is found that in this case the conversion 
factor has a different value which is independent of the beam size. In the limit for 
infinitely large beam sizes, the values of the conversion factors for both techniques are 
equal. 

The application of the theory discussed in this paper to actual observations will 
be discussed in a forthcoming paper. This will include a comparison between optical 
and radio observations. 

Key words: Methods: data analysis — ISM: general 



1 INTRODUCTION 

The accurate measurement of angular diameters is a long standing problem. This problem is pertinent to the study of planetary 
nebulae, Hll regions, galaxies and other exten ded s ources. Nevertheless, onl y few papers dedicated to t his p roblem can be 
found in the li terature: Mezger & H ender son (1967), Panagia & Walmsley (197J), Bedding & Zijlstra (1994), Schneider & 
Buckley ( [l996| ), and Wellman et al. ( |l997| ). 

Several methods are in general use to determine angular diameters. One method determines the angular diameter of 
the source, based on the Full Width at Half Maximum (fwhm) of a two-dimensional gaussian fitted to the observed surface 
brightness distribution (also called profile in this paper) in a least-squares sense. This method is usually called gaussian 
deconvolution and will be explained in more detail below. The second method is basically identical to the first, except that 
it determines the fwhm using the second moment of the profile instead of a gaussian fit. To discriminate it from the first 
method, it will be called second moment deconvolution. Both methods have the disadvantage that they yield a result that has 
no well-defined physical meaning. Hence, a conversion factor is needed to translate the result into something meaningful. In 
nebular research this usually is the Stromgren radius. This conversion factor depends on the method being used, the intrinsic 
surface brightness profile of the source and the resolution of the observation. 

The theory to determine the conversion factor is derived in this paper in several stages. First, in Section H some basic 
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assumptions and definitions will be given. Next, in Section W the full analytical theory for gaussian deconvolution will be 
given. In Section W this theory will be used to construct an alternative algorithm to determine the fwhm of an observed 
surface brightness distribution. In Section g the theory will be given for second moment deconvolution. In Section o the main 
conclusions will be presented. In Appendix N some additional proofs will be presented and finally in Appendix M the most 
important symbols that have been used will be defined. The appendices are only available in electronic form. 

2 ASSUMPTIONS AND DEFINITIONS 

The methods discussed in this article can be applied to observations at any wavelength. More in particular, they are valid for 
optical, infrared and radio observations. The resolution of these observations is usually characterized by the size of the beam 
profile for radio data or by the size of the point spread function for optical and infrared data. Throughout the paper the term 
'beam' will be used and it will be implicitly understood that it can also mean 'point spread function' where appropriate. It 
will be assumed that the beam can be approximated by a gaussian. In this paper the intrinsic surface brightness profile will 
be defined as the profile that would be observed with a perfect instrument (i.e. an instrument with infinite resolving power). 
For simplicity it will be assumed throughout this paper that both the surface brightness distribution of the nebula and the 
beam are circularly symmetric. This is a rather severe restriction; nebulae rarely are circular, and also for radio observations 
the beam usually is elliptical. However, this simplified case already yields interesting results which can be applied to actual 
data. 

As was already remarked, a conversion factor is needed to translate the fwhm diameter yielded by the gaussian or second 
moment deconvolution method into a Stromgren diameter. In this paper the Stromgren radius of the nebula will be denoted 
by rs. In the rest of this paper it will be assumed that the true diameter of the nebula is 0d = 2rs. The measured fwhm of 
the nebula will be denoted by "I> and the fwhm of the beam by $b- Throughout the paper the deconvolved fwhm diameter 
$d will be used, which is defined by 



$d = ^$2-$2. (1) 

This quantity is also commonly called the gaussian diameter. It should not be confused with the fwhm of the deconvolved 
or intrinsic profile, which in general will be difi'erent. The conversion factor to obtain the true angular diameter from the 
deconvolved fwhm can now be defined as 

ed=7*d ^ 7 = ed/$d = 2rs/$d. (2) 

This conversion factor is a function of the resolution of the observation, or to be more precise, of the ratio of the observed 
source diameter and the beam size. Hence an independent parameter (3 is chosen, which is defined as 

/3 = $d/<E>b. (3) 

In the following sections more details will be given of the techniques that have been used to calculate the conversion factors, 
both for the gaussian and the second moment deconvolution method. 

3 DETERMINING DIAMETERS USING GAUSSIAN FITS 

In this section the conversion factor will be studied for the case where the fwhm of the observed surface brightness profile is 
determined using gaussian fits. An implicit equation will be derived from which the value of 7(/3) can be solved for arbitrary 
p. This equation will also be used to derive the first terms of a Taylor series expansion of 7(/3) near /3 = 0. This implicit 
equation is derived in three stages. First it will be shown how to derive the width of a gaussian fitted to an arbitrary profile 
in Section 3J|. Next an expression for a surface brightness profile convolved with a gaussian of arbitrary size will be derived 



in Section 3.2. Finally in Section 3.3 these results will be combined to determine the conversion factor. In this process various 
lemma's will be used which can be found in Appendix |A[ 

Note 1: in the remainder of this paper all functions are implicitly assumed to be circularly symmetric. 

Note 2: in the remainder of this paper constraints have to be imposed on the surface brightness profile /(r). When c„ is 
defined as 

oo 

c„ = 27r/'/(r-)r"+Mr, (4) 



the following constraint can be formulated: 

(1) < fir) < F for all r £ [0, cxa), such that co exists and co > 0. 

This condition is sufficient for the theory derived in Section 3.1. For the theory in Section 3.2 (up to Eq. Gil) an additional 
constraint has to be imposed 

(2a) C2n exists for all n 6 No, with cj" = o{n) (n —^ oo). 
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It should be noted that the fact that Cn exists automatically implies that all Ci exist with i < n, as is shown in Lemma hi 



In the remainder of Section 3.2 



and all subsequent sections, this constraint needs to be restricted to 

(2b) C2n exists for all n G No, with c„'" — o{n^) {n — + co). 

It is also worth noting that all profiles which fulfill the following condition 

(2c) /(r) = for all r > Ro, 

(where Ro is some arbitrary outer radius) and which additionally fulfill condition (1), automatically fulfill condition (2b). 
This is proven in Lemma bl Since for all observational data conditions (1) and (2c) are automatically fulfilled, the theory 
presented in this paper is valid for all observed profiles. This neglects the fact that noise might cause some pixels to have 
negative flux values. In general this poses no problem however, as will be discussed in a forthcoming paper (van Hoof 
1999). One should also note that for a perfect gaussian, condition (2b) is not fulflUed. 

Note 3: in the remainder of this paper it is assumed that rs — 1. 

3.1 Fitting a gaussian to a surface brightness profile 

In this section the formulae to approximate a surface brightness profile /(r) with a two dimensional gaussian in a least-squares 
sense will be derived. The gaussian will be written as 

_ 2 

g{r;a,p) = ae ^^ with a > 0, p > 0. 

In order to determine a least squares fit to the surface brightness profile, the integral over the quadratic residuals x^ needs to 
be minimized. Hence the following two equations need to be solved: 

2-K oo 

^^' - ^ ' ''lf{r)-g{r;a,p)frdrd^ = 0, (5) 



dp dp , 

b 



and 

27r oo 

1^ = |- / [[fir)-g{r;a,p)frdrd^ = 0. (6) 

oa oa J J 



Eq. (H) will be evaluated first. After integration over (p and reversal of the order of differentiation and integration one gets 

oo 

-47r| [/(r) - 5(r; a,p) ] ^^^^^rdr = 0. 



One can easily prove that this is allowed using Lemma ^ and^, provided condition (1) is fulfilled. Since 

d9ir;a,p) _ _2-pr^ 
op 

one finds, after division by the constants, 
[ f{r) - ae"'"'' 1 r^e-^''\dr = 0, 



and thus 



f(r) r' e-P" dr = a r' e"^'"' dr = -^. (7) 

J 8p2 



For Eq. (|6|) one finds after a similar derivation that 

/(r)re-P'"'dr = ^. (8) 



Hence 

oc 

f{r)re~^^ dr — 2p f(r)r e~^^ dr, 



and thus the width of the gaussian fitted to the surface brightness profile can be found by solving Eq. (bI) . 
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Hp) = / .f{r) (1 - '2pr) r e'^" Ar = 0. (9) 



If the substitution s = pr^ is used, one can write 

CXD 

I{v) = ^j /(\/Vp ) (1 - 2s) e-'As = 0. (10) 



Once the width has been determined, the height of the gaussian can be solved from either Eq. (ph or Eq. (fef). Using Eq. (W 
one can easily see that 

lim7(p) = |^ > 0. 

P— >0 ATT 

Substituting a Taylor series expansion of / in Eq. (lid) one can prove that for p ^ oo the following is true 

/(p) = -/'"'(o)^^^±igi±^p-t-+o(p-t-i) < 0, 

where / (0) with n G No denotes the first non-zero derivative of /. Its value must be positive due to constraint (1). Because 
7 is a continuous function in p, its limiting behavior assures that at least one solution of Eq. (y_0|) must exist. Since 

^ = ^ >0, 
da'^ p 

this solution can either be a saddle point or a minimum depending on the sign of 

oo 
9V o « M/ 4 



27T- r^r-e-''" dr. (11) 

op p J dr 



From this equation it is clear that for monotonically decreasing /(r) all solutions of Eq. (hol) must be minima. Given the fact 
that x^ and its derivatives are continuous functions in both a and p, this implies that oidy one minimum can exist and thus 
the solution of Eq. (hol) must be unique. In general however, this is not the case. 



The constant surface brightness disk 

Eq. (y^ usually gives rise to implicit equations which can only be solved numerically. For a constant surface brightness disk 
the intrinsic surface brightness profile is very simple and a reasonably simple expression for the fwhm of the fitted gaussian 
and the conversion factor can de derived. 

I fo r < Ts 

/M= ' - (12) 

10 r > rs 

Without loss of generality it can be assumed that /o = I/tt and rg = 1. One can easily see that this profile satisfies constraints 
(1) and (2c). Substituting Eq. ( |l2[ ) in Eq. (hll) one can also prove that a u niq ue fit must exist. Hence, to determine the width 
of the gaussian fitted to the unconvolved profile, one may substitute Eq. (hd) in Eq. (hol) : 

p 

-^ //o(l-2s)e-Ms = ^ -!- [ (2p + 1) e-^ - 1] = 0. 
2p J Zpn ■' 



Multiplying this equation by 2p7r and taking the natural logarithm yields the following implicit equation for p 

p = ln(2p + 1) = 1.2564312 ... => ${p) = 1.4855024 ... => 7(00) = 1.3463458 . . . (13) 

3.2 An expression for the convolved surface brightness profile 

In this section an expression for an intrinsic surface brightness profile convolved with a gaussian of arbitrary size will be 
derived. This way a relation between the intrinsic and the observed profile will be found. The intrinsic (unconvolved) surface 
brightness profile will be written as /(r) and the convolved profile as fc{r). Hence 

00 00 00 00 

Mx,y)^ j j f{x',y')gix-x',y-y'-a,p)dx'dy' = a J j /(a;',y') e"''<"""'''e-^'""'''''da;'dy'. (14) 



The fact that the gaussian representing the beam should be normalized to 1 implies a — p/ir. In order to evaluate Eq. ( |14[ ) 
the exponential functions will be replaced by their Taylor series in the following way: 
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— p{x — x') —px —px' Ipxx' —px —px' \ ^ yZjpdjX j 



^'E 



(15) 



Since it can be proved easily that the series expansion in Eq. dla) is absolutely and uniformly convergent with an infinite 
radius of convergence both in x and x', it is allowed to substitute Eq. (|lfj ) in Eq. (|lj). Hence 



OO CXD 



Ux,y)^a / / f{x',y')e-^^ e'"^' e"''" e"'"'' 



E 



{2pxx' 



E 



Cipyy'Y 



Ax 'Ay' = 



OO OO 



-p{x~+y ) 



f{x',y')e 



(n — i)\ i\ 



n=0 ! = 



-Ax 'Ay' 



When X', yi are changed to polar coordinates, this equation can be rewritten as 

OO 27r 

fc[x,y) = ae-''<"'+"'' / / /(r') e^*""" V {2pY f^ ^f" 'f,'., r'" cos"" V sinV r' Aip' Ar' 



n=0 i=0 



(n-i) 



!i! 



(16) 



When conditions (1) and (2a) are fulfilled, the order of the integration and summation may be reversed (see Lemma |g) and 
one finds 



„n — 2„,z 



/c(a;,y)=ae-''<"'+^') / /(r') e-^"'" V (2p)"r'"+' V -^ ^ / cos"" V sinV dv?' dr'. (17) 

■'q n=0 i=0 ^ -' ^1 

The inner integral is well known and yields a non-zero result only when both n — i and i are even. In this case the result is 

27r 



2n-2^ . 2^ ■. „ (2n - 2^) ! (2^) ! 

cos (/9sm ydi/9 = 27r — -. 

2-'" n\ (n — i)! i\ 



When Eq. (|18|) is substituted in Eq. (hTf) and all the odd terms in n and i are omitted one finds 



fc{x,y) =, 



oi-^+v^) I ^^^/^„-P'■'^ V^ ^o„^2"„/2"+lV^ a:'" ''y'' ^_ (2n - 2i)! (2i)! ^^, ^ 



7/(-')e--''E(2p)"V^"-''E 



■ 27r 

(2n-2J)!(2i)! 2^" n! (n - i)! i! 



(18) 



1 " ^2n-2i 2i 

2n+i 1 v^ a; y 






27r(2e 



OO 2ri "■ / \ 



= 27rae 



/OO 2ri 
/(r')e--''5]^r'^"+i(xVy^)"dr'. 



When x, y are changed to polar coordinates and the order of integration and summation is reversed one gets (Lemma proves 
that this is aUowed) 



/c(r) = 27rae-P'''y ^r^" / f(r')r'^''+^e-'"''''Ar'. 
^-^ n'. I 

When the generalized n-th moment of f{r) is defined as 

OO 

Cn{p) = 2tt f(r) r"+"'^e~'"' dr, 



Eq. (uM can be rewritten to the following expression for the convolved surface brightness profile 

OO , , 2n 

, , _pr.2 ^-^ C2„ (p) P 2n 

/c(0=«e \^ — T . 



(19) 



(20) 



(21) 



In Lemma H it will be proven that this expression is absolutely and uniformly convergent for all p,r € [0, oo), provided that 
conditions (1) and (2a) are fulfilled. An alternative formulation for Eq. (Elf) can be derived which will be used further on. 
When the exponential in Eq. (GOl) is expanded in a Taylor series in p (which can easily be shown to have an infinite radius of 
convergence) one finds 
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6 P. A.M. van Hoof 

c4p) = 27r / /(r) r"+^ ^ (^-if^EI^Ar. (22) 



fc! 

"o '--=0 



When conditions (1) and (2b) are fulfilled, the order of integration and summation may be reversed (Lemma W), hence 

oo 
oo k r 



fc=0 







When Eq. (Ml) is used, one can write 

°° fc 

c4p) = E(-1)'|!C"+2^- (23) 

fc=0 

In Lemma H it will be proven that this series is absolutely and uniformly convergent for all p G [0, oo), provided condition (2b) 
is fulfilled. Hence one can substitute Eq. (E3|) in Eq. (Ell) 

°o 2n 2n oo k 

Mr)=ae--^ ^ ^ ^ (_!)'= |_ e,„+,,. (24) 

Lemma and Lemma H imply that this series is absolutely and uniformly convergent for all p G [0, oo) and r G [0, oo), provided 
conditions (1) and (2b) are fulfilled. Hence the summation may be reordered in any way. The following transformation will 
be used: n ~* k' , k ^ (n — k'), which yields (after dropping the primes) 

/e(r) = ae-- ^ C2„p" ^(-1)-'= -i^^. (25) 

n = fc=0 ^ '' ' 

Closer inspection reveals that the inner summation is closely related to the Laguerre polynomials L„. Hence one can write 

oo ^ 

Mr) = ae-^-' ^(_i)"£H!l^i„(p,2)^ (2g) 



3.3 The width of the convolved surface brightness profile 

Now an expression for the convolved profile has been found, the next task is to determine the width p' of a gaussian fitted to 
this profile. This is equivalent to measuring the fwhm of an observed profile. The result can be found by substituting Eq. ( |25[ ) 
in Eq. (|lO|). This is allowed since Eq. (bsl) is absolutely and uniformly convergent for all p G [0, oo) and r G [0, oo). The result 
is 

~ ~ \ k 

P^ ds = 0. 



OO oo 

^ f 1 /■ oo 71 

2^y(l-2.)e-Vc(V^)d. = ^ 2^yU-2.)e-ae-^E-"^"EMr-V-fc)lfc!HP 
When one abbreviates (" = -^ + 1 and multiplies with 2p/a, one can write 

oo 

/ ^ ^ /^ ■t\k 

(l_2.)e-^=E ^2„p"E(^l)"-'' fa_feVfc,2 ^''d. = 0. (27) 

Q n=0 fe = ^ '' ' 

Provided condition (2b) is fulfilled, the order of summation and integration may be reversed (see Lemma bt), hence 

CXD 

E C2„P"E('1)"" (n^Vlp /(1^2.).''e-^M. = 0. 



n=0 fe = 

It is well known that 

oo 

e dx = — —■ 

pTl + l 







I' 



Hence 



sn-fc (C-1)' 



{n-ky.kP 



k\ rfc + 1)! 



^fc + 1 (^fe4 



„_fe[C-2(fc + l)](C-l)' 



E«"pE(-i) f„,_MiM2 7iTr-2— ^^^ =0 ^ E=2-P E(-i) fn-Mifc. FiT^ = o 



{n-k)\k\ C* 



E^E(-ir-M;)[c-2(fc + i)](<^ = o. (28) 

n=0 ' fe=0 ^ ^ ' 
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Measuring angular diameters of extended sources - Theory 7 

Usually the width of a gaussian is given as the Full Width at Half Maximum (fwhm) value $. One can easily see that the 
FWHM for the beam $b is given by 

$b(p) = 2.Ji^. (29) 

The deconvolved fwhm $d of a fitted gaussian can be defined as 



*d(p',p) = ^-J'2(p')-*^(p). (30) 

Substituting Eq. (S) in Eq. ^) gives 



When one recalls the definition for /3 in Eq. (H) a relation between ^ and 13 can be obtained 



$d(p',p) ^ ln2,^ ^, / ln2 



When Eq. (bl) and Eq. (Bll) are combined, a relation between p and 7 can be found 



-<«^5^^/(C^^Vi^ - .(«^-,V)^=«. ,33, 

When Eqs. (||) and (||) are substituted in Eq. (||) one finds 
" C2„7'"(/3)/?'"ln"2^, ,,„_fc/n\ ,^2 ^,, (/32 + 1)" _^ 



^ C2„-J KP)P m . ^ ^_^^„_, ^ (^2 _ 2fc) 



fe=0 



fcy^'^ '(/32 + 2)*+2 



If it is assumed that /3 > this equation can be multiplied by (/3^ + 2)^//3^. This yields the following implicit equation for 

7(/3) 

E """^"If '"" ' E (-1)"-^ ( : ) (/^^ - 2fe) /^-- (g;^) ' = 0. (34) 

Now one can determine 7(0) by taking the limit /? — > of the left-hand side of Eq. (M). This yields 



co-C27'(0)ln2 = ^ ^(0) = ^-£^. (35) 

This yields the following expression for the deconvolved fwhm. 



c&d(0,0) = 2W^^i^. (36) 

V Co 

So it can be seen that the deconvolved fwhm of the source in the limit for very large beams is fully determined by two simple 
integrals over the intrinsic surface brightness profile. When Eq. ( pq ) is compared with Eq. (M) one can see that this limiting 
value simply is the fwhm derived from the second moment of the unconvolved profile. In Section H it will be shown that this 
is also the deconvolved fwhm derived from the second moment of a profile convolved with a beam of arbitrary width. 

Eq. ( |34[ ) can also be used to determine a value for 7(/3) for arbitrary (3. It can easily be solved numerically using a 
Newton- Raphson scheme. It can be shown that the Taylor series of 7(/3) contains only even powers of /3. Hence one may write 

oc 
1=0 

It should be noted that in general this expansion need not have an infinite radius of convergence, hence this expression is 
only valid sufficiently close to /3 = 0. Now the 2n-th power of this series can be calculated, which will have the same radius 
of convergence as the original series 

7^"(/3) = E»2rV' with sr'^So"; ffr = -^E[2M2n + l)-2n<;2.flg:k- (37) 

32 I 1 \ *: 



When Eq. (pTf is substituted in Eq. (M) one finds 



C2nin 2 ■r-^ {2n) V^ / , xn-fc / ^ \ , 2 r,i \ a2n+2l~2 ( P +1 



V-^ C2n -in A sr-^ (2n) V^ / 1 \n-fc / '' \ / ^2 r, 7 \ o2n+2i-2 P + ^ \ n ^oo^ 

E^^^EsL ^E(-l) U)(/5-2fc)/3 (1=0. (38) 

n = 1 = k = ^ ^ \l^ • / 
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towards 
observer 




Figure 1. The geometry for a constant volume emissivity shell is shown. Two lines of sight through the shell are also shown. The surface 
brightness for these lines of sight is proportional to the length of the thick parts of the lines. 

Taking the 2nd derivative, the 4th derivative, etc. . . to /3 of Eq. (BSf) and subsequently taking the limit /3 ^ 0, new equations 
are obtained from which pj: 54i etc. . . can be solved. Without proof the following result will be given 



7(/?) 



Co 



Cain 2 



2ci - C0C4 „2 , 2OC2 - I8C0C2C4 + 7Cocl - 2CoC2C6 „4 \ , ^,r,6 



Acl 



■r + 



2,24 



/3M+O(/?«)(/3^0). 



(39) 



Again it is noted that this series expansion in general need not have an infinite radius of convergence, hence it is only valid 
sufficiently close to /9 = 0. 



3.4 Application of the theory 

3.4.1 The constant surface brightness disk 

The coefficients C271 for the constant surface brightness disk can be calculated by substituting Eq. (113) in Eq. 



This gives 



C2n = 27r 



1 

/I 2n+l 



Ar — 



2n + 2 



Substituting this result in Eq. (B9|) yields (where additional terms are given without proof) 
"2~ /, 1 ^2 5 ^4 43 . 1367 o 51817 



liP) 



'-l^^i^' 



360'^ 12960^ 544320 



10 , 57119^i2|^Q(^i4)(^^Q)^ 



P'" + 



In 2 V 6" '36' 360" ' 12960" 544320" ' 653184 

Inspection of this expression suggests that the radius of convergence B is approximately B ~ 1 



(40) 



3.4.2 The constant volume emissivity shell 

The geometry for the constant volume emissivity shell is shown in Fig. hi Using this figure one can easily understand that the 
surface brightness for any given line of sight is proportional to the length of that part of the line that actually passes through 
the nebula. Hence for a constant emissivity shell the surface brightness profile is given by 

{2> ( Vr'i — r-2 — ^r? — r^j r<rs_<rs 
"iju^/r'i - r^ n<r<r^ (41) 

ri < Ts < r 

Without loss of generality it can be assumed that jv ~ j^ and r^ = 1. One can easily see that these profiles satisfy conditions (1) 
and (2c). First an analytic expression for C2„ will be deduced by substituting Eq. (Ull) in Eq. (M) 



C2n — 



f ^ ' i \ 

r " v 1 — r^dr — / r " \/rf — r^dr 



These integrals can be solved easily by using the substitution t = r^ and give 
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fc=0 



1-r 



2n + 3 



(42) 



k J 2n + 3-2k' 
The coefficients in the hmiting case for an infinitely thin shell can also be calculated. After normalization to co = 1 one finds 



C2r 



-EM 



,k+n / n 



2n + 3 



k I 2n + 3-2k' 



Substituting Eq. (Ua) in Eq. (bs) gives 



7(0) 



5(1 - rf) 
2(1 -rf) In 2' 



(43) 



(44) 



When Eq. (|43) is substituted in Eq. (B9) the following expression for the spherical case ri = is found (again giving additional 
terms without proof): 



-.(R\ / 5 , ^r2,5 4 1157 ,6 , 2226013 ,8 42074119 ,uA , ^,^i2w. n^ 
^(^)^V^h^l'"7^ +i5^ -11319^ +24720696^ " 519134616^ j + ^^^ )(^^°)- 

When Eq. MM is substituted in Eq. ( p9| ) one finds for the limiting case ri — *■ 1 

7(/?) = 



1 1«2^ 59 4 128 6 70117 o 20675393 10 \^^..i2w. „^ 
21n2V^^5^ +350^ " 875^ +539000^ ^ T^ET^HT^/^ ) + OW ) (/? - 0). 



175175000' 
Inspection of both expressions suggests that the radius of convergence B is approximately B ~ 1. 



(45) 



(46) 



4 AN ALTERNATIVE ALGORITHM FOR DETERMINING THE FWHM DIAMETER 

In this section, an expression will be derived that can be used to determine the fwhm diameter of an arbitrary surface 
brightness profile. This method will constitute an alternative algorithm to determine the fwhm diameter which is fully 
equivalent to a gaussian fit algorithm. The derivation will start with Eq. (p4): 



y^ ^^,, C2n7^"(/3)ln"2 " ^ ^^, ( n 



{0^ - 2k) P^ 



/3" + l 
/32 + 2 



0. 



(47) 



This expression will be transformed into an expression for the conversion factor of the intrinsic surface brightness profile 
by taking the limit /3 -^ oo. This will be done by substituting /3 = 1/e and subsequently taking the limit e — > 0. First the 
innermost summation (which will be called Sn) will be evaluated. 



s.-Ti-^f[l 



EM) 



i-2fe^.--^l + ^^ 



l + 2e^ 



EM 



''='5^) {l-2ke^)e-^" 



{I + e^)^{--l)' 2' e^' 



k / n 



(l-2fce^)e" 



l + 5^(-l)'2'-^^' 



^(-1)*= ? (l-2fce^).--^a,,e^ 



(48) 



The Taylor expansion of (1 + 2e^) ^ has a radius of convergence of |\/2, hence this derivation is only valid for < e < |\/2- 
The coefficients ak,i have been defined such that: 



l + ^(-l)'2'"'e^ 



a-k,i- 



(49) 



= '^a.k,ie^'- => afc.o = 1; at,; = - ^ [m(fc + 1) - i] (-1)"" 2'""^ 

i = m = l 

This expression for a^^j will yield a polynomial in k of degree /, which will be written as (see Lemma |lO| ): 

' fc' '"■^ k' 

ak,i =^di,,k' = {-!)''— + ^di,,k' = (-1)'— +fefe,i for allfc,iGNo and where 0° ee 1. 

1=0 ■ 1=0 

Here bk,i represents the terms in k with powers less than I (if any). It can easily be shown that the series given by Eq. (Ufl) is 
absolutely and uniformly convergent within its radius of convergence i \/2- Hence it is allowed to substitute the last expression 
in Eq. (Uq) and reorder the summation: 



^"-EM)M? (l-2fce^).--X 



{-^)%+bk. 



k 



= EM 

fc=0 
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When the limit e — > is taken, terms with I > n will tend to zero, hence the following is true: 






limSn^limy (-l)M 7 

fc = ^ 

When the order of the summation is reversed, this gives: 

i=o fe=o ^ ' ■ ;=o fc=o ^ ' 

n n y \ Ti— 1 n y \ 

+ E^"""E(-i)Mr)^^.' -2 E^"'"'''E(-i)'(r)fc''M. (50) 

i=0 fe=0 ^ ^ Z=0 *,— ^ ^ 

When the well known result 



E(-l)'° ( fc ) '^^ = (-1)"'^! S„,i for aU n,/ G No, / < n and where O" 
fc=0 ^ ^ 



is used, one can see that Eq. (BO) can be written as 

" 2i~2n "-1 o 2i-2n+2 n„| 

lim S„ = limE(-l)' -IT (-l)"'^' ^-' + E(-l)'^' ^^1^ (-l)"'^! ^".'+1 = 1 + T^T^ = 2n + 1. (51) 

The last two summations in Eq. ( JSOJ ) evaluate to zero because bk,i only contains powers of k up to fc'~^ or less. When Eq. ( |5l[ ) 
is substituted in Eq. (WTI) one finds the following result for the limit /3 — > cxa 

E(-l)"(2n + l) --"'^7""^ =0. (52) 

n = 

Which can be used to determine the limiting value of the conversion factor for any arbitrary surface brightness distribution. 
In the context of this section, an arbitrary surface brightness profile can mean an arbitrary observed profile, i.e. an intrinsic 
profile convolved with an arbitrary beam profile. To indicate this difference with previous definitions, primes have been added 
to the radial moments. When Eq. (H) is used, this expression can be rewritten into an expression which yields the fwhm of 
the observed profile directly: 

V^^ ixn2n + l , /41n2\" 

E(^l) ^^^2„' ^^1 =0. (53) 

n=0 ■ ^ ' 

This implicit equation uses only the radial moments of the observed profile and can be solved easily using a Newton-Raphson 
scheme. 



5 DETERMINING DIAMETERS USING SECOND MOMENTS 

It is well known that the fwhm of a gaussian is related to the second moment of a gaussian profile through: 



$,=2j£ii^. (54) 

V CO 

This formula is widely used to calculate the fwhm of an arbitrary profile. In general however, the result will not be identical 
to the FWHM derived from a gaussian fit. To distinguish the two values a subscript 2 has been used. Since the definition of 
the FWHM is not identical, also the value for the conversion factor will be different. One can define (using Eq. 091) 



2 



4>d = ^$i-<l>g(p) => 72(/3) = ^ = 



ln2(^-i 
Co' p 



(55) 



To distinguish between the values of the radial moments for the unconvolved and the convolved profile, a prime has been used 
in the latter case. Now an expression for the radial moments of the convolved profile will be derived. Substituting Eq. (Eq) in 
Eq. (0) one finds 

C2.' = 2.aEc2./E(-l)''"7fc:^ /-""'""^e--=dr. (56) 

fc=0 1=0 ^ 

Here the order of summation and integration has been reversed; in Lemma nil it is proven that this is allowed provided 
condition (2b) is fulfilled. The integral can evaluated easily using the substitution s = r^ and yields 
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{k-i)\m pi+^+^ p'^ ' ^-' ' (k-iy.ip' 

The fact that the gaussian representing the beam should be normahzed to 1 imphes a = p/ir. If an.i is defined such that 
^ a„,i I' = {l + n)...{l + 2){l + 1) = ^^ '^J^'' for all n,Z e No, and where 0° = 1 ^ 

i=0 

On.o = n\, a„,n = 1, ffln.i = na„_i,i + a„_i,i_i for all n G No, 1 < J < n — 1 
one can write 

fc=0 i=0 ^ ' 1=0 fc=0 i=0 i=0 ^ ^ 

The innermost summation is well known and yields non-zero results only when i > k. Since i < n, this also implies k < n. 
Hence one can write 



k 



C2: 



'-T.^E(-^)'-TrEi-^y {:u'^E<''7^ 



k\ '^^ ' \l ~ "^ "'"pn-k- 

i^k 1=0 ^ ' k=0 



The following result will be postulated for the coefficients b„^k 

This expression has been checked numerically as was found to be correct for all n,k < 19. This makes it very plausible that 
it is correct for all n and k. Substituting this result yields the following expression 

7 ;MM2 ~;r3T: =^ Co' = Co, C2' = I-C2, C4' = ^-H hC4, ••• (57) 

{n — k)\k\-' p"~k p p^ p 

Substituting Eq. ( |57| ) in Eq. (|5q), one finds 

'C2' 1 



72(/?) 



In 2 

Co' p 



ln2(i + ^-i 
p Co p 



7(0). (58) 



C2 In 2 



Thus it has been proven that the conversion factor for second moment deconvolution is independent of beam size and equal 
to the conversion factor for gaussian deconvolution in the limit for infinitely large beams. 

6 CONCLUSIONS AND FUTURE WORK 

In this work conversion factors have been determined to convert the deconvolved fwhm of a partially resolved nebula to its 
true diameter. This conversion factor depends on the fwhm of the beam and the intrinsic surface brightness distribution of 
the source. All work in this paper has been restricted to circularly symmetric surface brightness distributions and beams. The 
following results were obtained. 

(i) An implicit equation has been derived which can be used to determine the conversion factor given the intrinsic surface 
brightness distribution, the measured fwhm and the beam size. 

(ii) From this implicit equation, various explicit expressions have been derived, which give the conversion factor in cases 
where the beam size is larger than the source. 

(iii) Finally the implicit equation is used to construct an alternative algorithm for determining the fwhm of an arbitrary 
observed surface brightness distribution. 

(iv) The FWHM derived with gaussian deconvolution is in general not equal to the fwhm derived with second moment 
deconvolution. Hence the conversion factors will also be different for both methods. The use of second moment deconvolution 
is studied for the first time in this paper and it is found that the conversion factor is independent of the beam size in this 
case. In the limit for infinitely large beam sizes, the values of the conversion factors for both methods are equal. 

In a forthcoming publication the application of this theory to actual observations will be discussed (van Hoof 1999). Particular 
attention will be given to the limitations of the gaussian and second moment deconvolution method. Also a new method for 
deconvolving angular diameters will be presented. 
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APPENDIX A: ADDITIONAL PROOFS 

In this section additional proofs for the existence and convergence of certain integrals and series wiU be presented. 

Lemma 1 When c„ exists, then also all Ci exist with i < n. 
PROOF. One can write 

CX2 1 CXD CXD 

c, = 27t[ f(r) r'+'dr = 2-k f }{r) r'+'^dr + 2tt f f{r) r'+^dr < ^^ + 2'k f /(r) r"+^dr < 

1 1 

1 
<%F + Cn-2-K f(r) r"+'^dr < vrF + c„. 



Since f(r)r^^^ is positive everywhere. Hence the integral Ci exists for all i < n. 

Lemma 2 For f{r) for which conditions (1) and (2c) are fulfilled, condition (2b) is also fulfilled. 
PROOF. One can write 

oo Ro Ro 

C2„ = 2^ /" /(r) r2"+Mr ^ 27v f f(r) r^"+^dr < 2n F f r^^+^dr = ^^rI"+\ 



This proves that all C2n with n G No exist. Now one can write 
lim cli: < (JLI_RI-+A "^ ^ Urn (^-^] '^\l = Rl ^ lim cT = i?o = o{n^) 

n-too yn +1 J n-foo y n + 1 J n^oo 

and thus condition (2b) is fulfilled. 

Lemma 3 The integral In.k.i exists for all n, fc G No and all I (I'M provided condition (1) is fulfilled. The integral In,k,o exists 
for all n £ No and all fc £ N, provided that C2n exists and condition (1) is fulfilled. 

PROOF. To prove the first part one can simply write 

oo oo Ik 

Ir,,k,i= f''ir)g\r;a,p)r^"+'^dr<aF'' e''"'' r^"+Mr = ^ "'^ . 



To prove the second part one can write 

oo oc 

In.k.o = J .f'ir) r'"+Mr < F"'' J fir) r^"+'dr = ^ C2„, 



and hence the integral exists if C2n exists. 



Lemma 4 The integral I„,k,i is uniformly convergent for all a,p £ (0, oo), for any n,k £ No and any I £ N, provided 
condition (1) is fulfilled. 

PROOF. Lemma H proves that the integral In.k.i exists. Now one can use 

Y Y 

In,k,i = / f'ir)g\r; a,py"+'dr < o'f' J e-'^^\^"+'dr. 



X X 

If the substitution s = r^ is used, one finds 



a F f _i„„ „ , a F 



Since IpY^ > 0, one may write 
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a'F 



n + l Z-^ 



MpX 



2\n-i -IpX-' 



"■^•'-2(/p)"+i^(n-z)!' 

Since n\/{n - i)\ < nl and (IpX^)"'' < (IpX^)"- for X > {lp)~K one may write 

/„.M<5(;^);rTT(- + l)-!(/p^) e = ^^^ Xe . 

It is well known that Iimx-»cx) X^"e~''''^ = for any n G No, I G N and p G (0, cxa). Hence one can always find a finite solution 
^ ^ (^p)"' for the inequality 

I r' I ^ ^ \^2n -;px^ ^ 2/pe 

(n + Ij! a' _r '^ 

for all e > 0, n, fc G No, I G N and a,p £ (0, oo). This proves that In.k,i is uniformly convergent. 



Lemma 5 In the expression given in Eq. \ln) the order of summation and integration may be reversed for arbitrary values 
of p and r, provided that conditions (1) and (2a) are fulfilled. 

PROOF. Since the terms in Eq. (hd) have alternating signs, the following partial sum will be considered 



^-,« - E E 



^f[r')e-^^' [2pY 



„n — 1^ % 



{n ~ i)\ i\ 



r' cos ip' sin tp' r' 



di^'dr' 



After changing x, y to radial coordinates, one can write 



In,r<Y.Y. h (2P)"7;7^ /(^') '■'"^'d^'d-' = E E / 2- (2P)"^ " fir') r'"+Mr' 



EE2^« 



(2pr)" f n 



f{r')r'"+^Ar' < ^^2™ 



(2pr)" (n 



where for the last step condition (1) has been used. Since it is assumed that C2n exists for all n G No, this implies that c„ 
exists for all n G No (Lemma hi). Hence In,r exists for arbitrary A'' and R. This also proves that \miR-^oc In,r exists. For 
arbitrary R the following is true 



(2pr)" 



(2pr)" 



■ c„ 

n = ■"' i = ^ ' n = '" n=0 



(4pr)" 



W<^2™^^c„^ r^ =^2™^^c„2" = ^2™^:i^c„ = 5]a„ 



In order for limjv^oo In,r to exist, the following must be true 



lim a,/" — lim 

n — >cxD n — -oo 



(4pr)" 
iyra- -^ c„ 



i/r. 



< lim 



■ (4per)" 
27ra-^^ —c„ 



l/n 



i/„4per i/„ 



— lim (27ra) / " "^" ' ^ v " _ jjj^j 

n — ►cxD 72 n — ►cxD 7X 



^^C^" < 1. 



Here the well known result n! > e~"n" for all n > has been used. In order for the last inequality to be true for arbitrary 
values of p and r, the constraint Cn = o{n) (n -^ cxd) has to be imposed. This proves that limjv^oo In,r exists. Since the 
argument given above is also valid in the limit R -^ oo, this also proves that limjv,i?^oo In,r exists. This completes the proof 
that the order of integration and summation in Eq. (hq) may be reversed. 



Lemma 6 The series given in Eq. ]2l\) is absolutely and uniformly convergent for all r G [0, cxa) and p G [0,00), provided 
that conditions (1) and (2a) are fulfilled. 

PROOF. First it will be proven that these conditions imply [cn(p)]^ = o{n) (n -^ cxj) for all p G [0,00). Since /(r) > 
everywhere, one can write 



C2n(p) = 27r / f{r)r 



2n + l -pr^ 



dr < 2tt f{r) r "'^^dr = C2n. 


This proves that C2n{p) exists for all n G No and all p G [0, 00 ). It also implies 

[c„(p)]^''" < cl/" = o(n) (n -^ 00), 



for all p G [0,cx3), which completes the first step. The summation in Eq. (Ell) is a Taylor series, which will have an infinite 
radius of convergence when the following condition is fulfilled 
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lim 



C2n{p)p^ 



l/{2n) 



< lim 



C2n(p)(pe) 



l/(2n) 



lim El[c,„^p)]m2n) 
n — -oo n 



Here the well known result n\ > e^"n" for all n > has been used. For the last condition to hold for any arbitrary value 
of p, it must be true that [c2n{p)Y ~ o{n) (n — > cxd) or [c„(p)]"'^'" = o{n) (n — > cxi), which has been proven to be true. 

Any Taylor series is absolutely and uniformly convergent within its radius of convergence. Since the gaussian in Eq. ( pi|) is 
bounded for all p,r £ [0,00), this also proves that the expression in Eq. (Gil) is absolutely and uniformly convergent for all 
p,r e [0, 00). 



Lemma 7 In the expression given in Eq. IkS) the order of summation and integration may be reversed for arbitrary values 
of p, provided that conditions (1) and (2b) are fulfilled. 



PROOF. Since the terms in Eq. (E2|) have alternating signs, the following partial sum will be considered 



= E 



o / 1\l<: P J^r \ n+2fe + l 

27!- (-1) -j^finr 



dr = ^ 27r — / /(r) r^'^'^'^^'^dr <^'^-^ -^ Cn+2k = ^ Ok, 

k=0 ' n fe=0 ' *.— 



where for the leist inequality condition (1) has been used. Since it is assumed that C2n exists for all n G No, this implies that 
Cn exists for all n G No (Lemma hi). From this it is clear that Ik,r exists for arbitrary values of K and R. It is also clear that 
lim_R_>oo Ik,r exists for any value of K. In order for limx^oo Ik,r to exist, the following must be true 



lim aj/*" = lim ( 27r ■— Cn+2k 



i/fc 



< lim 

k — *oo 



27r ^^. C„+2k 



l/k 



lim(2.)V^^,V.^^. 



lim -rrc^'k < 1- 

fc— >oo fc 



Here the well known result fc! > e *fc'° for all fe > has been used. In order for the last inequality to be true for arbitrary 
values of p, the constraint c^). — o{k) {k -^ 00) has to be imposed. Using the substitution n — 2k, this can be written as 
c„ = o{n2) (n -^ cxj). This completes the proof that liniA'^oo Ik.r exists. Since the argument given above is also valid in 
the limit R —> 00 this also proves that limK.R^oo Ik.r exists. This completes the proof that the order of integration and 
summation in Eq. (G2|) may be reversed. 



Lemma 8 The series given in Eq. Isjj is absolutely and uniformly convergent for all p £ [0,oo), provided condition (2b) is 
fulfilled. 



PROOF. The summation in Eq. (|23| ) is a Taylor series in p, which will have an infinite radius of convergence when the 
following criterion is fulfilled 






l/k 



V 6 1/fc e l/k 



Here the well known result fc! > e~^k^ for all fc > has been used. Condition (2b) assures that this is true, which can be 
seen by making the transformation n = 2fe. This proves that the series has an infinite radius of convergence. Since any Taylor 
series is absolutely and uniformly convergent within its radius of convergence, this completes the proof. 

Lemma 9 In the expression given in Eq. (127} the order of summation and integration may be reversed for arbitrary values 
of p, provided that condition (2b) is fulfilled. 

PROOF. Since the terms in Eq. (ETI) don't always have the same sign, the following partial sum will be considered 



E 



(l-2s)e-«=C2„p"^(-l)"- 



(n-fc)!fc!2 



ds < ^ C2n p" ^ 



(C-i)" 

(n-fc)!fc!2 



(2s + 1) s*"' e'^^'ds < 



JV 



<E^-p"E 



n=0 

JV 



(n-fc)!fc!2 



o (fc + l)! , k\ 



fc + 1 



C + 2fc + 2 

C2 ^^ (n-fc)!fc! 

n=0 ^ fc=0 ^ ' 



E^E 



, Y^ C2nP"(C + 2n + 2) ^ (n\_^ ca,. (2p)" (C + 2n + 2) _ ^ 



n=0 



C^ 



c^ 



Here the fact has been used that C > 1 (see Eq. K2|) and hence ((^ — 1)'° < (^''. Since it is assumed that C2n exists for all 71 G No, 
it is clear that Jjv.s exists for arbitrary values of A'' and S. It is also clear that lims^cx) In,s exists. In order for limjv_,oo In,s 
to exist, the following must be true 



© 1999 P.A.M. van Hoof 



Measuring angular diameters of extended sources - Theory 15 



lim a„'" = lim 

n — 'CxD n — *oo 



— lim 

71 — -OO 



C2„ (2p)" (C + 2n + 2) 



C2, 



n! C 

(C + 2n + 2) 



l/n 



< lim 



C2n(2pe)"(C + 2n + 2) 



C^ 



1/t7 



C^ 



1/n 



2pe i/n2pe 
= lim Co' < 1. 



Here the well known result n\ > e "n" for all n > has been used. In order for the last inequality to be true for arbitrary 
values of p, the constraint Cj^ = o{n) (n -^ oo) has to be imposed. Using the substitution 2n -^ n, this can be written 
as Cn = o{n^) (n — » oo). This completes the proof that limjv->oo In,s exists. Since the argument given above is also valid 
in the limit S ^^ oo this also proves that lim]v,s^oo In,s exists. This completes the proof that the order of integration and 
summation in Eq. ( P7| ) may be reversed. 

Lemma 10 The expression for Ok.i given in Eq. (udj yields a polynomial in k of degree I for all Z G No. 

PROOF. From the definition given in Eq. (Ml) it is clear that the lemma is correct for i = 0. If it is assumed that the 
lemma is correct for all I' < I then the last part of Eq. ([43) assures that a^^i is equal to the sum of I polynomials of degree 
I' — I — 1, 1 — 2, . . ., multiplied by a linear term in k. This yields a polynomial of degree max(i' + 1) = Z. Since the lemma is 
true for / = 0, this proves by induction that the lemma is true for all / G No. 

Now an expression for the highest order term of this polynomial will be derived. It is clear that the highest order 
contribution to ak.i only comes from the term with 7n = 1 in Eq. Wm). If one defines 



ak 



1,1 = ^ di,ik' 



one can deduce the following expression for dij from Eq. ( |49| 

-di-2.i-2 = 






n(-i)- 

Li=0 



do,o = (-l)'^ 



i=0 



Here bk,i is the lower order remainder of Ok.i which is a polynomial of degree Z — 1 or less in k. It can easily be checked that 
these expressions are also valid for A; = and Z = if one postulates 0° = 1 and 6^,0 = 0. 



Lemma 11 In the expression given in Eq. (pn) the order of summation and integration may be reversed for arbitrary values 
of p, provided that condition (2b) is fulfilled. 



PROOF. Since the terms in Eq. (M) have alternating signs, the following partial sum will be considered 



lk,r 



E 



27rac2fcp'°E ("1)' 



{k-l)\l\ 



2l + 2n+l -pr'' 

-r e 



dr = E 27ra C2k p^ E 



(k-l)\l^? 



2l + 2n + l -pr' 

re 



dr< 



<'^ TVa C2k p^ E 



(Z + n) 



E 



!=0 
k 



{k-iy.ip p'+ 



^I+r = E^-=2fcP E 

— n ^ 



{i + ny. 



(k~iy.ip 

1=0 ^ ' fe=0 



-E 



C2kP 



E 



k\ 



{l + n)\ 



(k-iy.U l\k\ 



< > C2fcP" 
fc=0 



^ , k\ {l + nr ^^ 



fc-n \k + ny 



k\ 



E HE 



C2fcP 



fc_„2'=(fc + n)^ 



k\ 



K 

E«fc- 



In this derivation it has been used that the gaussian used to convolve the profile should be normalized to 1, this implies 
a = p/tt. Since it is assumed that C2k exists for all Zs G No, this proves that Ik,r exists for arbitrary values of K and R. It 
also proves that limH_*oo Ik,r exists for arbitrary values of K. In order for IvaiK^x Ik,r to exist, the following must be true 



lim ttj, = lim 

k — >oo k — >oo 



C2fcP 



k-r.i'^ik + nY 



— lim c 



i/fe 

■2k 



k + n 



k\ 

n/k 



1/fe 



— lim 

fc — »oo 



C2k 



k + n 



(2p)* 



Zc! 



i/fe 



< lim 

k — >oo 



C2A; 



k + nY (2pe)* 



Zjfc 



l/k 



2pe 



fc— >oo ^ P / k fc — 'OO 

-fcy^fc 



1/fc 2pe 
'2fe k 



Here the well known result Zc! > e k for all Zc > has been used. In order for the last inequality to be true for arbitrary 
values of p, the constraint Cj^ — o{k) {k -^ oo) has to be imposed. Using the substitution n — 2k, this can be written as 



Cn ~ o{n^) (n —> oo). This completes the proof that lim/v'^oo Ik.r exists. Since the argument given above is also valid in 
the limit R ^ oo this also proves that \imK,R^oo Ik,r exists. This completes the proof that the order of integration and 
summation in Eq. (pw may be reversed. 
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APPENDIX B: SYMBOLS 

The following definitions for the symbols have been used. 

a -Measure for the height of the gaussian. 

Cn -Radial moments of the surface brightness profile. Defined in Eq. M). 

Cnip) -Generalized radial moments of the surface brightness profile. Defined in Eq. 

f(r) -The surface brightness profile of the nebula. 

F - Upper bound for the surface brightness profile. 

(?(r; a, p)- Two-dimensional gaussian profile. 

jv -The emissivity per unit frequency. 

Ln -Laguerre polynomial of degree n. 

p - Measure for the width of a gaussian. Related to the fwhm through Eq. 

r -Polar coordinate on the sky. Subscripts have the following meaning: 

i - the inner radius of the nebula, 

s - the Stromgren or outer radius of the nebula. 
x,y -Cartesian coordinates on the sky. 

/3 -Ratio of the deconvolved fwhm diameter of the nebula to the fwhm of the beam. 

7 -Factor to convert the deconvolved fwhm to the true angular diameter. Subscripts have the following meaning 

<none> - gaussian fits were used to determine the fwhm 

2 - second moments were used to determine the fwhm (Eq. p5[ ). 

/ - indicates the fitting function used to approximate 7(/3). 
r -The gamma function. 

Si J -Kronecker delta symbol. 

^ -Auxiliary variable, related to (3 through Eq. (621). 

Qd -True angular diameter of the nebula (Od = 2rs). 

ip -Polar coordinate on the sky. 

$ - The FWHM diameter of a profile. Subscripts have the following meaning: 

<none> - indicates the fwhm of a gaussian profile, 

b - indicates the fwhm of the beam, 

d - indicates the deconvolved fwhm of a profile (defined in Eq. BOl or Eq. JSq) , 

2 - indicates the fwhm of a profile obtained using second moments (Eq. p4^ . 
Other symbols have varying meanings as defined in the pertinent sections. 
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